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rS ■ Abstract 

We investigate the Balitsky-Kovchegov (BK) equation ioi D = 3 space-time dimensions, 
corresponding to one transverse coordinate, and we show that it can be solved analytically. 
The explicit solutions are found in the linear approximation and for the paricular cases 
when they do not depend on the impact parameter of the dipole or on its rapidity. It is 
shown, that in a general case different solutions are related by an infinite parameter group of 
transformations. The key observation is that the equation has the Kowalewskaya-Painleve 
property, which gives the possibility of reducing the non-linear problem to the solution of a 
linear integral equation. 
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1 Introduction 

In the framework of perturbative QCD the most general basis for the theoretical description 
of small-x processes is given by the BFKL approach [1], based on the gluon Reggeization. 
Taking into account also analyticity and unitarity constraints for the S'-matrix one can cal- 
culate the high energy asymptotics of scattering amplitudes with a good accuracy [2]. In 
particular, using this approach the kernel of the BFKL equation was obtained in the next-to- 
leading order (NLO) [3]. To avoid a strong renormalization scale dependence of its solution 



in the MS'-scheme one can chose the Brodsky-Lepage-Mackenzie (BLM) optimal scale set- 
ting within physical renormalization schemes [4]. It is important that the BLM procedure 
allows to preserve, approximately, some good properties of the leading order approximation, 
in particular, its conformal invariance [5]. Incorporation of the corrections related to the 
renormalization group symmetry [6] should make the predictions of the BFKL equation for 
the small-x physics more reliable. 

However, at asymptotically large energies the BFKL approach should be generalized, 
since it leads to a power-like growth of cross sections. The problem of the Froissart bound 
violation cannot be solved by the calculation of radiative corrections at any fixed order 
and requires other methods. One of the possibilities is to take into account the multi- 
reggeon contributions in the framework of the BKP equations [7] . A general framework for 
considering unitarization effects is the reggeon field theory, describing the interactions of 
reggeized gluons. A possible way of deriving such a field theory is based on the high-energy 
effective QCD action [8], generalizing the BFKL approach for multi-particle production and 
multi-Reggeon interactions. 

In a first step, the unitarization program can be realized by non-linear generalizations 
of the BFKL equation, related to the idea of saturation of parton densities [9]. The most 
important ingredient is the transition vertex for 2 — i> 4 reggeized gluons which, in the BFKL 
approach, has been calculated in [10] and further investigated in [11, 12]. Using this transition 
vertex, as a first step beyond the BFKL equation, one can formulate a nonlinear evolution 
equation which sums the QCD fan diagrams consisting of BFKL Green's functions [13]. In 
a future step towards finding solutions to the reggeon field theory, closed loops of BFKL 
Pomerons should be included and calculated. 

At present it is fashionable to describe the saturation in DIS within the color dipole 
approach [14], which in the target rest frame has an attractive and clear physical interpre- 
tation: the incoming 7* splits into a qq colour dipole long before the dipole interacts with 



the target. In this approach the cross section is presented as 



a^.{x,Q^) = j d\ j dz\^j.{r,z,Q^)\^adp{r,x), (1) 

where \1/^* (r, z, Q'^) is the photon wave function, z the quark momentum fraction, and r = |r| 
the transverse size of the colour dipole, r*= r*i — r*2; the transverse vectors r*i and r*2 are the 
quark and antiquark coordinates, and adp{r,x) is the dipole cross section. 



adpir,x) =2jd'h{l- 5(fi,f2;r)) . (2) 

Here 6 = (ri + r2)/2 is the impact parameter of the quark- antiquark system, Y = log(l/x), 
and S'(r*i,r2; Y) = (tr {U^{fi)U{f2)) /N^y is the S'-matrix for the elastic scattering , which 
can be presented as an average of the product of two Wilson lines 
U^{x_[_) = P exp[ig J dx~A^{x~ ,xj_)t°'] . 

There exist two different approaches for calculating the Y dependence of S{fi,f2;Y). 
The first one starts from an infinite hierarchy of coupled equations [15] for products of any 
given number of Wilson lines. Taking the limit of large Nc and considering the case of the 
target being a large nucleus, the equation for the two-line correlator decouples and takes a 
simple form [16]: 

dS{fi,f2,Y) asNc /• 2^ (ri-f2)2 ^ . ,^ ^ . 

- d r— --— —-[S{ri,r,Y)S{r,r2,Y) 



dY 27c^ J {fi — r)2(r2 — r)^ 

-S{n,f2,Y)]. (3) 

At small N = 1 — S, one can neglect terms nonlinear in A^, and this equation turns into the 
colour-dipole version of the BFKL equation for the dipole cross section A^. In [13] it has been 
shown that Eq.(3) can be considered as a special limit of the QCD fan diagram equation 
mentioned before: apart from taking the large- A'c limit, the dipole cross section has to be 
in the Mobius class of functions. An alternative approach for describing saturation is the 
Color Glass Condensate [17]. In this framework, the F-evolution is based upon functional 
equations. It is affirmed [18] that both approaches give identical results for the calculation 
of observables. 

Eq.(3) has attracted much interest. First, numerical solutions have been found for a 
6-independent approximation [19], and later on also for the full 6-dependent equation [20]. 
Several attempts have recently been made to find analytical solutions, at least in some 
approximation [21]. However, so far no exact analytical solution is known. Therefore it will 
be interesting to find a model in which the BK equation could be solved. 



In this paper we generalize the BK-equation (3) to an arbitrary dimension D of the 
space-time. For the case D = 2 + 1 - which corresponds to one transverse dimension - we 
shall present analytical solutions (the BFKL dynamics was investigated at D = 2 + 1 in 
the papers [22]). The key observation in finding these solutions is that the equation has 
the Kowalewskaya-Painleve property. We proceed in several steps of approximations. After 
writing down the equation in one transverse dimension (Section 2) we first discuss the linear 
approximation (Section 3). Then we turn to the nonlinear equation, beginning with the h- 
independent case (Section 4), where we find explicit solutions for arbitrary initial conditions. 
Passing on to the 6-dependence we first find all "stationary" (rapidity independent) solutions 
(Section 5), generalizing then to the rapidity dependent case (Section 6). We conclude with 
the investigation of the infinite parameter symmetry of the equation (Section 7) and its 
general solution (Section 8). A few details are put into the Appendix. 



2 The fan equation for 2+1 QCD 

Let us begin by writing down the fan equation for the S'-matrix describing the dipole evo- 
lution in the (2 + l)-dimensional QCD. The kernel of the BK equation (3) is equal to the 
probability density in the space of transverse coordinates and rapidity for the soft gluon 
emission by the quark-antiquark pair: 

asN, (ri - f^f 2a,iV, | / ^^i(e^^>-^T - e^^>-^^) | ' (4) 



27r2 (fi - f)2(f2 - f)2 ' ' ' y (27r)^-2 p 

at D = 4. It is clear that if the space-time dimension takes the value D = 4 + 2e, the BK 

equation preserves its form, but since 

(i^+^'^/c exp {ikr)k i r(l + e) / f \ 

- (5) 



(6) 



(27r)2+2^ k2 2tt n' {{7^)^+^ J ' 

the kernel changes according to 

asN, (fi - f2)2 asNc /r(l + e) \ V (n - f) (fs - f) 



2vr2 (fi - f)2(f2 - f)2 27r2 \ 71' J \^(fi-f)2(l+^) (f2-f)2(l+^)y 

Putting e = — 1/2 we obtain: 

'"^' - dpo {Sp2Po{y)Sp„p,{y) - Sp,p,{y)) , y = ^y , (7) 



dy J^^ -ru,-^p,poy.j--poPi..j -P2P.,.n.. 27r 

where the inequality p2 > Pi is imposed without any restriction of generality. Note that the 
scattering amplitude Np,^p^{ii) is related to the S'-matrix Sp^p^{y) by the identity 

Np,pM = ^-Sp,pM- (8) 



We introduce the dipole impact parameter b and its size p 

= — ^ — , p = p2- pi, p = po- Pi ■ (9) 

Then the fan equation can be presented as follows 

dS{b, p, y) fP ^^, / ^^ ^p'^_ ^^ ^^ ^^^ _ p__l^ ^^ ^^ _ ^^^^ ^^ \ ^ ^^^^ 



dy k ' V 2'^ r,..-v^ 2 

or, in a simpler way, 

(^ + P) ^(^' P' ^) = // ^P''5(^ - ^' P - P'' ^) ^(^ + S^' P'' ^) • (11) 

Now we go to the mixed representation (6, p) 

-^e''^.(6,p,l/) (12) 

-joo ZVTi 

with the Laplace transform 

/•oo 

#,p,t/) = / rfpe-''^5(6,p,t/). (13) 

JO 

In this representation p is the differential operator 

d 
P sib, p,y) = -— s{b, p, y) = -ds{b, p, y) . (14) 

The convolution of two functions in the p-representation 

Fuip)= rdp'F,ip-p')F,ip') (15) 

Jo 

is reduced to the product of their Mellin transforms 



fuip) = flip) f2{p) , fip) = / dpe-^^Fip) . (16) 

Jo 

Moreover, a more general relation holds 

fuip) = (d^fiip)) {d'f2{p)) (17) 

for the Mellin transform of the function 

Fi'^ip) = ['dp' {-p')\-p + pyF,{p - p') F,ip') . (18) 

Jo 

Let us introduce the anti-normal ordering A^ of operators depending on the two momenta Pr 

(r = 1,2) and the corresponding derivatives dr = d/d{pr), according to the rules 

N{prds) = N{dsPr) = dsPr = Prdg + Ks ■ (19) 



Then one can verify that the fan equation in the mixed representation (6, p) can be written 
in the form 

d~~d~) ^*^^'^'^) ^ ii%'^ (s{b+-d2,pi,y)s{b- -di,p2,y)] 1, (20) 

where it is imphed that, after performing the transformation of the differential operator to 
the anti-normal form, one omits all terms with the differential operators due to the relations 

d[d^l = (}, k,l>l 
and, at the end, puts pi = P2 = P- 

Note that one can consider a modified equation for the S'-matrix in the case of a non-zero 
temperature T in the t-channel (see [23]). In this case S{p,y) is assumed to be a periodic 
function of pk 

Sp2Pi{y) = 'S'p2Pi(i/)iPfc^Pfc+i/T • (21) 



3 The BFKL equation in 2 + 1 dimensional QCD 

To begin with, we consider the linearized case corresponding to the BFKL equation in 2 + 1 
dimensional space-time. Using the relation S{b, p, y) = 1—N{b, p, y) and neglecting nonlinear 
contributions in (11), one can obtain 

[h ^ ^ ^^^' '"' ^^ ^ r '^p' (^^^ " 2 "" ~ p'' y^ + ^(^ + ^^p- p'^ y)) ■ (22) 

This equation is simplified significantly in the case when N does not depend on b 

(^ + p\ N{p. y) = ^l' dp'N{p', y) . (23) 

For the Mellin transform 



oo 



n(p, y) = I dp e-P^ N{p, y) , N{p, y) 9{p) = H ^ e^ n{p, y) (24) 

JO J-ioo iTXl 



we obtain the differential equation 



d d 2 

;^-^)^(P'^) = -^(P'^)- (25) 



Its solution with the initial condition n(p, 0) = no(p) is 



{y + pfnQ{y + p) 
n{p, y) = 2 • (26) 



The second order pole at p = provides a linear dependence of N{p, y) at large p 

N{p,y)c:^py^no{y), (27) 

which leads to a violation of the S'-matrix unitarity aX y ^ for any initial condition. Note 
that Eq.(23) has a stationary (y-independent) solution N{p,y) = cp, with arbitrary c. But 
this solution is not acceptable from the physical point of view. 

Let us turn now to Eq.(22) for the amplitude depending on b. Performing the Mellin and 
Fourier transformations in the variables p and b, respectively, 

N{b,p,y)9ip) = r ^.e^ T ^ e^«^n(g,p,|/) , 



oo /-oo 



n{q,p,y) = dpe-^^ dbe-'"^ N{b, p,y),) (28) 

Jo J-oo 

we obtain 

Its general solution is 

/ X f{p + y;q) /gn^ 

n{q,P,y)= ^^,^^, , (30) 

with the arbitrary function f{p + y',q)- For the initial condition n{q,p,0) = no{q,p) the 
result can be written as follows 



{4:{y + p f + q'^)no {q, y + p) 
4p2 _|_ qi 



n{q,P,y)= '^''^'> ::>T'' " ^ (31) 



where 



Therefore 



no{q,p)= dpe""" dbe'''"' N{b, p,0) . (32) 

Jo J-oo 



dp ^^ r dq j,i, (4(y + pf + q^) 

2 I /,2 



N{b,p,y) = l ^e^ I ^e'^' 

-ioo 2m J-oo 27r 4p^ + q 



X dp'e-p'^P+yn db'e''^'''N{b',p',0). (33) 

JO J-oo 

By introducing the Green function Gy{b, p;b', p',0) according to the definition 

fOO fOO 

N{b, p,y)= dp' e-^ y / db' Gy{b, p; b', p') N{b', p', 0) , (34) 



'0 

we obtain for it 

Gyib,p;b',p') =5{b-b')5{p-p') 

+0{p-p) TT^ M2|/cos + sm , (35) 

J-oo Ztt \ Z q 2 / 



so that the Green function for the BFKL Pomeron is a polynomial in y 

Gy{b,p;b',p') =6{b-b')6{p-p') 

+oiP - p') E {y^i^ -^' + ^^) + V^(& -y + ^^)) • (36) 

Therefore the y-evolution of the initial dipole distribution Np^p^ (0) in the (/)2, pi)-representation 
is simple: 

N,,M = e-''''Np,pM + y r '^P^ (^""'"^p^A.lo) + e-yp--Np,pM) 

Jpl 
+y' r '^P'^ /"' dpi^~^'''~'''^'Np>^p,^ (0). (37) 



'pi •>pl 

For large y ^ 1/p, where in (30) one can neglect p in comparison with y, the solution in 
the coordinate representation is 

N{b, p, y) = r ^ e-^ !i!lM^( V + q') no{q, y) . (38) 

J-oo ZTT Iq 

In particular, for large b^ pwe obtain the above discussed expression N{b^ p, y) — py"^ ^o(y)- 
On the other hand, at b <^ p one has 

7V(M.!;)-r|^5EM?) 2,. ,„„(„), (39) 

J-oo ZTT g 



Note, that Eq. (29) admits stationary solutions 



m 

4p2 _|_ g2 



^(9^p^y)=7i^^ (40) 



with an arbitrary function f{q). For N{b,p,y) we obtain correspondingly 

N{b,p,y) =a(6 + |)-a(6-|) (41) 

with an arbitrary amplitude a. Actually, the existence of such solutions can be seen from Eq. 
(22). They are given by a difference of the amplitudes depending on transverse coordinates of 
quark pi and antiquark p2, respectively. Splitting of the dipole does not lead to any evolution 
of such a solution: new quarks and antiquarks on the edges of the dipole have the same 
transverse coordinates, and therefore their contributions cancel in the linear approximation. 



4 Solutions depending only on the dipole size 

Next let us investigate the simple case when the solution does not depend on b, which, 
formally, corresponds to the asymptotics of s{b,p,y) at 6 ^ oo 

lim s{b,p,y) = s{p,y). (42) 

0— »oo 

This function describes the S'-matrix for the scattering of the dipole with the size p on the 
hadron system with a large radius R^ p, provided that the impact parameter b is restricted 
to the interval 

R:>b::^ p . 

In this case the fan equation is significantly simplified. In terms of variables 

^=2^y-p), V = ij^{y + P) (43) 

it takes the form: 

^5(e,^) = s^(e,^). (44) 

In accordance with the fact that this is a degenerated Riccati equation, this equation has 
the Kowalewskaya-Painleve property [24], and its solution is a simple meromorphic function 

K^.,)^^. (45) 

Returning to the p-representation, we obtain the following result 

-loo ZTii i^s[p + y)) — y 
where s{p) is the Laplace transform of the S'-matrix at the initial moment y = 

s{p)= / dpe-^S{p,Q). (47) 

Assuming the colour transparency property in the form 

hm S{p, 0) = 1 - cp^ + ... , 7 > , c> , (48) 

one can calculate the function s{p,y) at large y or p 

Therefore, if < 7 < 1, we obtain in the region p ~ y^^"^ <^ y the scaling behaviour 

s{p, y) ^ — — Tv^ — r-[zz , (50) 

p + cT{l + 7)1/^ ^ 



which leads to the S'-matrix 

S{p,y)-e-''^'^^+^^y"''p (51) 

in the region y^^^'p ~ 1. It means, that at large y the distribution in p is proportional to a 
smeared (5-function S{p). 

If, on the other hand, 7 > 1 (for example, when limp_,o S{p, 0) = 1 — cp^) the S'-matrix 
tends to unity in the large region p <^ y^~^, which corresponds to the vanishing of the dipole 
interaction. In the particular case when 7 = 1 and 

lira S{p,0) = I- cp + ..., oO, (52) 

the asymptotic behaviour of the S'-matrix at large y is especially simple: 

\im S{p,y)=e-^^. (53) 

It corresponds to the fact that in the 2 -|- 1-dimensional QCD the fan equation has the 
stationary solutions 

S{p, y) = e-'P , (54) 

corresponding to the conservation of the total length of dipoles. This expression is similar to 
the Maxwell distribution in energies for molecules in a thermal equilibrium. In our case the 
energy is proportional to the total length p of dipoles. Thus, there is an analogy between 
the Boltzman and BK equations. 



5 Stationary solutions of the fan equation 

The stationary fan equation 

0= / dP0Wp2P0^P0p\ ~ ^p2Pl) (^^) 

J pi 
has the following solution 

5 _ S{p2) , . 

where S{p) is an arbitrary function. In a linear approximation it corresponds to the ampli- 
tude N{h, p, y) (41) for the BFKL equation, being the sum of two functions depending on pi 
and p2. In order to show that Sp^p^ is a general solution, we expand the equation in terms 
of small fluctuations ASp^p-^ around the solution 

^^P2P1 ~ ^P2pl ~ ^P2P11 ['^') 

10 



and we obtain, in the linear approximation, 

0= / dpo ISp^pf^ ASp^^pj^ + ASp^pf^ Sp^^pj^ — ASp2p-^j . (58) 

By introducing the relative correction 5Sp^p^ 

^Sp2Pl = ^P2pl^^p2pl 5 (59) 

we simplify the above equation 

rp2 , — —- ,^ . 

0=/ dp^ibSp^p^^bSp^p^-bSp^pA . (60) 

J p\ ^ ' 

There is an obvious solution of this equation 

^Sp2Px = s(p2) - s(pi), (61) 

where s{p) is an arbitrary function. This correction, however, represents only a redefinition 
of the function S{p) in the expression for Sp^p-^^. In fact, by introducing the new variables 
p = P2 — Pi, b = (p2 + Pi)/2 we can reduce the above equation for bSp,^p^ to the stationary 
equation for N{b,p,y), which has only the solutions bSp^p^, as it was shown in Section 3. 

The stationary fan equation has the Kowalewskaya-Painleve property. Namely, if we the 
solution in the series 

oo 

Sp2P, = T.^n{h)p''; ao = l, (62) 

n=0 

the first coefficient ai{h) is an arbitrary function of p with the following asymptotic behaviour 
at 6 ^ oo 

«i(&) = i:|- (63) 

r=l " 

Indeed, the above solution of the stationary equation can be written as follows 

oo /!+ p Y^ 

fe=l V 2(b+dk) I 

where d]~ and r^ are arbitrary complex parameters. Therefore, the above function ai(6) is 

In contrast to this, as a consequence of the Kowalewskaya-Penleve property, for a given ai(6) 
the solution Sp^p^ is obtained in an unique way. For example, if we start from the scaling 
ansatz 

aM = \. (66) 

11 



the stationary solution is 



'b + ^Y 



,2 ^2 / „ „3\ ^3 



S.»HF^I -l+«^^^+^ + ^^ + 0(/). (67) 



6-f / 6 2 62 ' \i2 ' 6 / 63 



We can also write solutions of the stationary equation which do not have any singularities 
at real pi,2- For example, 



6 Scaling solution of the non-stationary equation 

Let us now consider the general case when s{b,p,y) depends on 6, and let us introduce the 
new variables 

V = ^{y + p), e = ^{r]) --{y-p), (69) 

where ^(r^) is an arbitrary function which later on will be derived from the requirement that 
s{b,p,y) satisfies initial conditions at y = 0. In terms of the new variables the fan equation 

looks as follows 

d — 

— ^—S(b,e,ri)= hm A^ 

(*(" + Kll + <' + ^'<*>' ^) • "• '"> ^<'' - 3 (^ + <' + ^'<"'" s;) • -• *)) ' ■ 

(70) 

Its solution can be searched in the form of an expansion in powers of 1/e: 

1 °° ~ ~ c (b n) 

S{b, e,v) = - + Y: KS{b, e, V) , KS{b, e, V) = 4# • (71) 

f r=l ^ 

Here Ci{b,ri) turns out to be an arbitrary function, fixed by initial conditions for s{b,p,y) at 
y = 0. According to the fan equation the residues Cr{b, rj), near 6 = oo, are analytic functions 
with the following asymptotic behaviour: 

CrM = lt'-^. (72) 

In particular, for the first correction AiS{b,e,ri) to the pole 1/e we obtain the differential 
equation 

(-^^-'^)A,Sib,e,v) = 0. (73) 

Its solution is 

Ai5(6,6,r^) = ^i^, (74) 

12 



where 

ci(^^) = lE-I^ (75) 

" s=0 '^ 

is an arbitrary function which is analytic and vanishing at 6 ^ oo. It should be chosen to 
satisfy the initial conditions for s{h,p,y) at y = 0. 

Therefore, in the expansion of S'(6, e, Vj) in powers of 1/e, there is a phenomenon similar 
to the Kowalewskaya-Painleve property for the integrable differential equations. Namely, 
one can construct the general solution in the class of meromorphic functions by fixing the 
functions ^{rj) and Ci{b,ri) from the initial conditions for s{b,p,y). The next correction 
A2S{b,e,ri) satisfies the equation 

(_|_!)a,s(.,,„) = 

(76) 
where, in the last term, it is implied that after the differentiation of Ci the removed operator 
b + ig^ is substituted by the expression 1 +^'{r]2) in the corresponding place. Therefore we 
obtain 

A2S{b,e,v) = '^, (77) 

where 

C2(^r^) = iVci(6+- — ,7^0-Ci(&-^^,r72)i„.,--^Ci(6+- — ,r^)(l+e'(r72))|^^^^. 

(78) 
In an analogous way one can calculate the residues Cr{b, rj) in terms of ^{tj) and Ck{b, rf) with 
k < r. 

Let us now consider the solution of the fan equation in the scaling regime 

e < 1 , 6 > 1 , e6 ~ 1, (79) 

for fixed tj. In this case the equation (70) is significantly simplified 

-|- S{b, t) = hm N (s{b + M^, ei) S{b - u^, e^)] 1 , (80) 

where the parameter u = u{ri) is defined by the expression 

u = \{l+e{v)), (81) 

13 



and the solution has the following expansion at large eb 

^(&,e) = i + ^ + .... (82) 

Here the function v = v{ri) is fixed by the initial conditions for s{b, p, y). The first two terms 
on the rhs of (82) determine all the other coefficients of the expansion 

S(b, e) = - y^ , ," ; an = 1, ai = v. (83) 

The reccurence relations for the coefficients a„ are given in the Appendix, and one obtains 

It is easy to see that the substitution 

^(&,e)^i-s(6,^), e-^2up (85) 

turns (80) into (20) for the stationary (y-independent) case. Therefore we have the remark- 
able result: the general solution of the non- stationary equation at e6 ~ 1 can be expressed 
in terms of the stationary solution depending on p/b. Using the stationary solution obtained 
in the previous section (see (67)) one can derive 

V 

1 poo / h ~\- ^ \ 2tt 

S(^,)^-l dpe-^^) . (86) 



7 Infinite parameter symmetry of the fan equation 

In this section we consider another ansatz: 

^P.P. (y) = s,,p. (y) = e^^^^^^-^^p^^^y ^ , (87) 

where R{p) and T{p) are arbitrary functions. Putting this ansatz into Eq.(7) it is easy to 
verify that it is indeed a non- stationary solution, provided that 

For small fluctuations 5Np^p^ (y) around this solution 

Sp,pAy)^Sp,pMil + SNp,p,iy)) (89) 

14 



we obtain, inserting (89) into (7) and linearizing the result in 5Np,^p^{y)) 
dSNp,pM _ r,„ RiPo) 



dy Jp, "^^'n^ i^Np^Poiy) + SNp,p,iy) - 6Np^pM) . (90) 

We have "zero mode" solutions of this equation 

5Np,p,iy) = 5R{p2) - 5T{p^) + y {5x{p2) - 6x{pi)) , (91) 

where SR{p) and ST{p) are arbitrary functions related to Sx{p) as follows 

Sxip) = r dpo^ (SRipo) - 6Tipo)) . (92) 

Jo T(po) 

Introducing the new integration variable T(p) 

T(p)^/.p„f|^, (93) 

we obtain the following linear equation for 5Np^p^{y) = Ny^y-^{y): 

^^h^M = Q dTo (iVT.Tod/) + Nr.rAy) - Nr.rAv)) , (94) 

which coincides with the BFKL equation having the general solution given by (33) with the 

substitution 

T -I- T 

2 - i 1 , 6 ^ ; T,. = T{pr) . 

Let us factorize Sp^p^ (y) as follows 

Sp2Pi (y) - Sp^p, (y) Sp^p^ (y) . (95) 

Then we obtain the following equation for Sp^p^{y)) 

9sp2P,{y) fp^ R{po) I I . I . I .. /n^^ 

-' dpo-——[Sp^p^[y)sp^p^[y)-Sp,^p^[y)). (96) 



dy J PI T{po) 

Therefore the fan equation has a renormalization group property. Namely, Sp^p-^^ (y) satis- 
fies the same equation as Sp^p^ly), but in the new variables T, defined by (93). As a result, 
we obtain the following "Backlund" transformation between different solutions of the fan 
equation 

e(---)^5i:L(l/)= |^e(-^-"^^^4'k(l/) (97) 

-'12IP1J 

for arbitrary R{p) and T{p), provided that the new variable satisfies 

T(p,^/'..fM. ,98, 
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The function T(p) is chosen to grow with p and to have the additional constraint 

T(±oo) = ±oo. (99) 

We can once more perform the above transformation: 

,(T.-T.).42)^^(^) = |H^e(^-^^)^5gK(y) , r. = r(T.) , (100) 

-'231 J-lJ 



where 



and verify its group property 



^ i?23(To) 



r(T)= / rfTo-^^, (101) 

-^231 ioj 



3^^^"^^^^^Si(^)= 144 e^'^''^^' 41(1/), r, = r(T(p,)). (102) 



Here 

^13(P2) Rl2{p2) R23i^2) 



(103) 



Tisipi) Tuipi) T23(Ti) 

rm,))^ f"\ir. MM. f,,,MeAMIMl. ,io4, 

J T23(To) ^ Ti2(po) T23(T(po)) 

In particular, for 

r(T(p)) = p (105) 

one obtains 

RUJP) i?23(T(p)) _ 

Ti2(p) T23(T(p)) • ^ ^ 

Generally, 

-Rl3(P2) ^ ^12 (P2) ^23(^2) , . . 

Tl3(pi) Ti2(pi) T23(T0 ^ ' ^ ^ 

and we thus obtain a family of solutions 

for arbitrary /(p). The symmetry group of the fan equation is a product of this abelian 
transformation of the dipole S'-matrix Sp^p^{y) and of the reparametrization group (98). 

Without any loss of generality it is convenient to write the action of the second group on 
S'p2pi {y) in the form 

^{p.-P.)ySi^) iy) = a{p^)a{p,)e^''-^^^y" Sr.rAy) . T, = T(p,) = r dpa\p) , (109) 
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where a{p) = J Ruip) /Tuip) ■ We can express this equation directly in terms of the trans- 
formation p' = T(p) 

e(P.-P0.5(i)^(^) = ^l^^^^e^^^'^^^yS^lrM ■ (110) 

It means that under the general covariance transformation p -^ T(p)the quantity 

T,,,M=e'^''^''^'S,,,M (111) 

transforms as a tensor with the spinor indices P2, Pi- The functional expressing the S'-matrix 
Sp2pi{y) through initial conditions 

^P.Pi(y) = 0(5p.Pi(O)) (112) 

should be invariant under the full symmetry group (together with the above multiplication 
of 'S'^jpi ^y ^^^ ratio f {P2) / f (pi)) ■ In particular we can use the above transformations to 
change in a convenient way the initial conditions 

SS.(0)=f|^sST.(0)^ (113) 

8 General solution of the fan equation 

The fan equation for the function 

Gp,pM=e^''-''^'Sp,pM (114) 

can be written as follows 

'"^' - dpoGp,,p^Xy)Gp,p,iy). (115) 



dy J PI 

The function Gp^p-^ly) can be considered as a matrix element of the infinite-dimensional 
matrix G{u) having a triangular form, in accordance with the fact that it is nonzero only 
for p2 > pi 

Gp,pM = {P2\G\pi) . (116) 

We can define the product of two matrices G°'{u) and G^{u) as the matrix G^{u) with the 
matrix element 

g''p,pM= rdpoG;,pMG\.pM- (117) 
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Correspondingly, the inverse matrix G ^{u) is defined as the matrix satisfying the relations 

S{p2-Pi)= dpoGp^p^Xy)iG )p^^p^{y)= dpo{G )p;po{y)Gp,,p,{y). (118) 

J pi J pi 

Its matrixelements (G^^) p.^p-^{y) are non-zero only for p2 > pi. Thus, the operators G{y) 
belong to the (solvable) group of triangular matrices. The matrix element of the unit element 
of this group is 6{p2 — pi). 

In the operator form the fan equation looks as follows 

'-^ = CAy). (119) 

and its general solution has the meromorphic property of Kowalewskaya-Painleve [24] 

Note that in accordance with the Ricatti theory the operator G^^ satisfies the free Newton 
equaton 

^G-i(y) = 0. (121) 

The corresponding solution for the S'-matrix Sp,^p-^ (y) written in the operator form is 

where p is the operator with the eigenvalues pi for the eigenfunctions \pi). This solution is 
similar to the general solution of a linear equation expressed in terms of its Green function. 
In a matrix form it can be expanded in a series as follows 

e ^ ^ ^P2pi\y) = '^P2pi(0) + 

/■P2 rp2 rPo' 

y dpoSp2po{0)Sp,p,{0) + y / dpo' dpoSp,p^^,{^) Sp,,po{0)Sp,p,{0) + ... . (123) 

jpi Jpi J pi 

In particular, according to above discussions one can simplify this general result as follows 



Sp.Pi iy) = e(>^(-)-(-))^ i^ , xip) = y; dpo (^i^ - 1 j , (124) 

providing that Sp^p^ (0) has the factorized form 

^P.Pi(0) = ^. (125) 



It is also obvious that the solutions Sp^p^ {y) for different Sp^p-^ (0) are related each to other by 
the above established group of transformations depending on two arbitrary functions f{p) 
and T(p) 



In particular, we can consider the case where the initial condition S't2Ti(0) is only a 
function of T2 — Ti, which gives us a possibility to construct the solution depending on 
three arbitrary functions 

^.P, (l/)-y_^ 2vr^ {s{p))-^-y' ^^^^^ 

where the parameter a is chosen from the condition that all singularities of the integrand 
are situated to the left of the integration contour in p and 



Jo 2 ' 2 



The solution S**--^'^'*-* can be considered as a particular case of the general solution transformed 
to a diagonal form by an appropriate similarity transformation. The function s{p) has the 
interpretation of the eigenvalue of 5'p2pi(0), and c^p^{p) is its eigenfunction. The similarity 
transformation is unitary (c*^^-* = c*^^)+ = l/c*-^-*) provided that /(p) is a phase: /* = f~^. 

In general the functions c^p\p) have a more complicated (non-exponential) dependence 
on p which enumerates the eigenvalues. This dependence can be found from the solution of 
the equation 

/•oo 

d Pi Sp2p, (0)cp, (p) = s{p) Cp^ (p) . (129) 



For the triangular matrix Sp^p^{()) the eigenfunction Cp{p) can be calculated easily, pro- 
vided that we substitute the continuous coordinate p by a finite number of points pk- In 
general, the behaviour of Sp^p^{y) at large y depends on the initial conditions encoded in the 
asymptotics of the functions Cp{p) at large p. Due to the colour transparancy property we 
have 

Sp,p, (0) < 1 , Jim^ Sp,p, (0) = 1 - cpl, 7 > , (130) 

where the coefficient c > is a function of h. At large y and fixed pr it is natural to expect 
an universal dependence of Sp^p^{y) from y. In particular, for large c the S-matrix should 
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have the form of a smeared 5(p2i)-function. It is related to the fact that in this case, as a 
result of the y-evolution, the average size of the dipole tends to zero, and one can neglect 
the fe-dependence of S (up to a common factor). In the pre-asymptotics a scaling behaviour 
of Sp^p^iy) is possible (cf. section 4). 

It is important that the solution of the Balitsky-Kovchegov equation satisfies also the 
linear integral equation 

Gp2Pi (y) = Gp^p, {0)+y dp Gp,^p^ (0) Gp^p, {y) . (131) 

•I pi 

In particular, the theory of the Fredholm equations can be applied to this equation. In this 
case, the perturbation theory in y should work rather well. For large y one can use a quasi- 
classical approximation. We hope to return to the investigation of the this linear equation 
in future publications. 



9 Conclusions 

In this paper we have obtained analytic solutions of a simplified version of the nonlinear BK 
equation, namely the BK equation in one transverse dimension. Our explicit formulae allow 
to study, for any initial condition, the rapidity evolution, both as a function of the dipole 
size and of the dipole position in impact parameter. 

We also have found several remarkable properties of the nonlinear equation: the general 
solution has the meromorphic Kowalewskaya-Painleve property, and the equation has sym- 
metries which are reminiscent of general covariance. Some of these results might be helpful 
in addressing also the more realistic case of two transverse dimensions. 
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Appendix 

Here we derive the reccurence relations for the coefficients a„ in (83). Inserting the 
expansion (83) in (80) one can obtain the relation 

Using the expansion 



oo 






T{k + r)T{l + r) 



'dx' -^ / r=0 " V I -^y-^ V'V \\ "-^Z -^ / r=0 

where 



r(r + i)r(A;)r(/) ' 

one can easily verify that oq = 1, and a\ is arbitrary. For the coefficients a„ with n > 1 we 
obtain the recurrence relation 

(n - l)a„ = 

ra— 1 n m— 1 n—ra 

Y: a, cij' (u--' + {-uT-') + E E ^'^ ^-'^ ^'^""^ E (-1)' cU-.+i cr-Tii . 

fc=l m.=2 A:=l /=0 

Evidently, the terms with odd powers of u are cancelled. 
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